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Abstract. We prove here essentially sharp L q (L r ) linear and bilinear esti- 
mates for the wave equations on Minkowski space where we assume the initial 
data possesses additional regularity with respect to fractional powers of the 
angular momentum operators Qij := x l dj — x^di. In this setting, the range of 
exponents (q, r) vastly improves over what is available for the wave equations 
based on translation invariant derivatives of the initial data, or uniform decay 
of the solution. 



1. Introduction: Classical Strichartz estimates, improvements for 
spherically symmetric data, and knapp counterexamples 

The aim of this work is to prove mixed Lebesgue space estimates for solutions 
to the linear wave equation on Minkowski space in a setting where the initial data 
is assumed to possess extra regularity with respect to weighted derivatives in the 
angular variable. These types of estimates arise naturally in applications to the 
scale invariant global existence theory of non-linear wave equations which do not 
possess certain "null" structures in their non-linearities. For example, in a compan- 
ion to this article, we use the estimates proved here to show global existence and 
scattering for the (4+1) Yang-Mills equations in the Lorentz gauge for a certain 
class of small, scale invariant initial data. 

All of the estimates we prove will be of "Strichartz type" , i.e. L q (L r ) space-time 
estimates for solutions of D<fi = 0. Due to the presence of extra weighted angular 
regularity, we will get a significant gain over the usual estimates for the wave equa- 
tions which are based solely on translation invariant derivatives of the initial data. 
What separates our estimates from the "classical" Strichartz estimates, is that they 
are not solely based on the uniform decay of solutions to the wave equation. In- 
stead, we will exploit a certain "wave packet" structure these solutions exhibit in 
radial coordinates. This allows us to decompose our waves into a sum of pulses, 
each of which remains coherent for all time. 1 These pulses can then be treated on 
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^This fact makes the type of wave packets we use here more simple to handle that those 
used in say 1181 . In fact, all of our estimates are proved by directly integrating in time, without 
using any complicated induction procedures or further analytical tools such as the TT* argument. 
However, there is a certain lack of coherence our wave packets exhibit, which seems to limit their 
usefulness to situations where one can use some angular regularity. This is because, strangely 
enough, our wave packets become more coherent for large times. For relatively small times, close 
to the "focusing" time of our wave packets, much of this coherence seems to get lost. See Remark 
14. 21 after the statement of Proposition ^, l| 
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an essentially individual basis. This works well to prove a certain weak 2 endpoint 
estimate for the range we are considering. We then interpolate our endpoint with 
the endpoint from JB] to obtain the full sharp (up to an e) set of estimates. 

We shall also consider multilinear type estimates which involve weighted angular 
regularity on one or any number of the factors. For this set of estimates, we again 
obtain a vast improvement over the "classical" multilinear estimates for the wave 
equation (e.g. compared with 0). We shall state and prove these estimates in 
the context of T. Tao's dual scale machine for generating multilinear estimates (see 
|16p. All of the estimates we prove here are sharp, up to an arbitrarily small loss 
of angular regularity, when tested against Knapp counterexamples. We will explain 
this in more detail in the sequel. 

After the proof of the estimates we describe here was discovered by the first 
author, a shorter proof was found by Igor Rodnianski for the case of 3 ^ n spatial 
dimensions. We have included this in an appendix to the present work and have 
elected to retain the discussion of our original proof in the main body of the pa- 
per because it includes the development of machinery that is interesting in its own 
right and is perhaps more flexible. Furthermore, our proof gives a lot of detailed 
and interesting information as to how wave propagation works on high spherical 
harmonics. In particular, we provide what seems to be a sharp localization of band 
limited Hankel transforms, the type of which has been studied by previous authors 
(see jH]). It is likely that the general procedure we employ here which uses this type 
phase space localization is applicable to other dispersive phenomena (e.g. wave or 
Schrodinger equations) on spherically symmetric backgrounds. 

We now begin with a brief discussion of the usual Strichartz estimates for the 
wave equation. The standard reference for this material at this point is the paper 
0. Let u\ be a unit frequency solution 3 to the wave equation on Minkowski space. 
By this we mean a function u\ such that: 

(1) □«! = o , S u P p{^(o),a^i(o)} c {iei 1 1 < iei < 2} . 

Without loss of generality, we may assume that u\ is of the form ui(t) = e u ^~ A fi, 
for some unit frequency function /1 of the spatial variable only. Then two of 
the most basic quantities which determine the space-time behavior of u\ are the 
following estimates: 

(2) ||e ltv ^/i|| L 2 - WflUl, (Energy Estimate) 

(3) || e 4tv/=S /i || i» < t 1 ^- || /1 || L i . (Dispersive Estimate) 

By interpolating between (J2J and ©, and using some standard duality arguments 
(specifically the TT* method), it is possible to show the following set of space-time 
estimates for u\: 



2 Not in the sense of real interpolation. By weak we simply mean an estimate which contains 
logarithmic divergences. 

3 For definitions of the various objects we present here, including Fourier transforms and mixed 
Lebesgue spaces see the next section. 

For more detail, see the notation list in Section 121 
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Theorem 1.1 ("Classical" Strichartz estimates including endpoints (see 0, 
and Let 3 n be the number of spatial dimensions, and let a = then 

the following estimate holds for 2 ^ q: 

(4) || e^Mi,^) < H/ilU, , 

where ^ + f- ^5 §, with the exception of the forbidden L 2 (L°°) endpoint on M 3+1 . 

A key facet of the proof of (J2J, is that it does not rely on any other property of 
the evolution operator e ** v/_A than the estimates @ and ©■ In other words, the 
proof only uses the conservation of energy and the uniform decay estimate Q. 
Furthermore, the estimate l@J is sharp in that one cannot improve the range of 
L q (L r ) indices stated above without replacing the L 2 norm on the right hand side 
of with something else. This can be seen as follows: Let us consider initial data 
sets fl = XB'i where \b> is the indicator function of a radially directed block of 
dimensions 1 x e x . . . x e, lying along the £i axis between 1/2 < £i < 2. Then a 
quick calculation using the integral formula: 

(5) e'^ftix) = Je 2 ™^+^XB*(0dH, 

shows that one has | e lt V-A ye ^| ^ e n-i Qn ^ ne S p ace -time region S% x : 

t = 0(e- 2 ) , 
t + xx = 0(1) , 
x' = Oie' 1 ) . 

Based on this one sees that: 

(6) e*<*-i-*>||A e |U" ~ e^WxstJw) < W^^fth^) - 

where a — Therefore, the condition 5 + 7 ^ f of Theorem 11.11 must be 

satisfied. We note here that initial data sets f{ are commonly referred to as Knapp 
counterexamples. 

With the above consideration in mind, it is natural to wonder if that somehow 
the initial data /1 were forced to be more evenly spread out along the various radial 
directions in Fourier space, then one could gain some improvement on the range 
of indices in However, any such improvement must somehow involve another 
mechanism than just the estimates @ and (J3J. This can be seen simply from the 
fact that both of these estimates are sharp even for spherically symmetric data (e.g. 
looking at waves which focus at lager and larger times). Even so, it has been known 
for some time that, with the help of additional arguments based on the specific form 
of the integral representation for e ltv/ ~ A , one can obtain a significant improvement 
over Q for spherically symmetric initial data (see for example |B] and [H]). To 
understand how this can happen, consider a unit frequency radially symmetric 
initial data set /°, such that f® is a smooth bump function of the radial variable 
|£| . Then by using an integration by parts and stationary phase argument on the 
integral representation Q5J, in conjunction with the phenomena of finite speed of 
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propagation it is not difficult to see that one has the asymptotics (because we can 
let derivatives fall on /°, of course this works for smooth enough Fourier data even 
if it is not spherically symmetric): 

Cm 



(7) 



|e^- A A°(r)| < 



|t|^(l+||t|-r|) 



M 



That is, at time t, e lt ^ A ff is essentially the indicator function of an O(l) spherical 
shell of radius t multiplied by the amplitude Based on this, one can easily 

computes that: 

II iiV-A fO II < \ 

Therefore, in order for us to have that the L 2 (L r ) norm of e lt ^~ A is finite, we 
need that i < (n — 1)(^ — ^), or equivalently, that 2( -" ^ 



n-2 



< r. This is a vast 



improvement over the requirement of ^ r coming from the Knapp coun- 

terexamples. The key to this improvement is that along with the uniform decay 
rate Q, the waves g^V-A jo are highly localized in physical space along the radial 
variable. 



Of course, an arbitrary spherically symmetric wave will not have the localiza- 
tion l|7|. However, it is possible to in a straightforward manner chop up a unit 
frequency spherically symmetric wave into a sum of pieces, each of which satisfy a 
time translated version of the asymptotic J2J. This is accomplished via a suitable 
physical space localization of the radially symmetric Fourier transform as follows: 
We begin by rewriting the integral formula (jSJ) for radially symmetric initial data 
/i as: 

2vr 



(8) 



7i(r) 



J^l(2nrp)fi(p)p^ dp 



where J n -i (y) is the Bessel function of order 



We now us the well known 



asymptotics for Bessel functions of relatively small order (see 117)1: 
(9) J^i(y) = < 



cos(y- 2^7r) -mi(y) 



sin(y - ^tt) • m 2 (y) 



for 1 ^ y , 

j^-m 3 (y), for ^ y < 1 . 

Here the function is C°°, and the remaining m, have asymptotic expansions: 



kV 



-2k 



m 2 (y) = 53 C; 



2,kV 



-2k-l 



^The asymptotic JJJ in the interior of the light-cone can be proved using integration by parts- 
stationary phase. In the exterior of the light cone, this asymptotic is true thanks to finite speed of 
propagation. Strictly speaking, the functions {f® (0), i\/—Af (0)} are not compactly supported. 
However, they are exponentially localized around the origin, so one may recover J7J in the exterior 
via weighted energy estimates. We will give an independent proof of this fact in a moment (see 
1141 ) which does not rely on weighted energy estimates. 
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as y — > oo. In other words, the functions mi(2irrp) and m2{2 / nrp) are C°° with 
derivatives in p uniformly bounded for all \ ^ p ^ 2 and 2 ^ r. Substituting the 
asymptotic (0 into the integral formula, we may assume without loss of generality 
that we are trying to bound integrals of the form: 

(10) I±(t,r) = — - / e 2 ^ t±r ^m ± (r,p) X{1/iA) (p)f 1 (p)dp, 

(1 + r)^~ J-oo 

where m is a smooth function with derivatives in p uniformly bounded for all 
^ r, and X(i/4,4) is a smooth bump function on the interval (t,4). It is now 
apparent that the integrals in l|l(J|) are essentially time translated inverse Fourier 
transforms of a one dimensional unit frequency function. Therefore, we can localize 
these integrals in physical space (i.e. the t±r variable) on an O(l) scale. This can 
be accomplished with the help of the so called (^-transform (see |2!), which is just a 
smoothed out redundant (over-sampled) version of the classical Shannon sampling. 
Since the function /i is compactly supported in the interval (0, 4), we may take its 
Fourier series development: 

(11) h{p) = c k e^ k " , p € (0,4) . 

k 

An important thing to notice here is that we can recover the L 2 of fi as a function 
on M. n in terms of the {ck}' 

(12) \\fi\\% ~ £ |c fc | 2 . 

k 

This can be seen from the Plancherel theorem and the fact that fi is unit frequency, 
so the volume part of the integral in Fourier space which comes from integrating 
over spheres is O(l). Sticking the series (|llf) into the the integrals l|10|) yields: 

(13) /*(*,»•) = ]T - g^- ^ftr) , 

^ (1 + r) — 



k 



where: 



^(t,r) = / e 2 ^ t±r +i^m ± (r,p) X{1/4A) (p)dp 



Integrating by parts as many times as necessary in the above formula, we see that 
we have the asymptotic: 

(14) |^(*,r)| < r - w 



l + \t±r + \\) M 



Using the expansion (|13|) and the asymptotic l|14|l we can directly compute that for 
2 < p: 



< 



< 



y * — ^) P (i +r) (»-i)(i-f) dr 

4- (l + |t±r+||)3 J 1 
\ck\ 



y / — w (i + r )<«-i)(i-f) dr , 

4- loo (1+ |i+f|-r)^ ' 
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The manipulation to get the last line above follows from Holders inequality and the 
fact that 2 p. By integrating each expression in this line term by term and using 
the inclusion £ 2 C £ p for 2 ^ p we arrive at the bound: 



||l S < E 



< 



„ (l+lf+fl) p(n - 1)( '-*\ 
E 



\c k ? 



w (i + |t+j|)>c-i)(i-i); 

Testing this last expression for L 2 in time we see that: 

(l+it+jr^HW)* 

-oo 4 



Now as long as 1 < 2(n — 1)(| — i), or equivalently 2 ^_^ < Vi we have the result: 

ll J± lli?(LS) £ E i c *i 

Using the characterization 1)12(1 . we have shown that: 



n-2 

2 



Proposition 1.2 (Unit frequency Strichartz estimates for spherically symmetric 
data). Let u\ be a unit frequency, spherically symmetric solution to the equation 
□«i = in 3 ^ n (spatial) dimensions, then the following space-time estimates 
hold: 

(15) \\uih HLn < |hi(0)|| i2 , 

where r satisfies the bound 2 ^ 1 _ 2 1 - ) < r. 



Interpolating (|15|l with the energy estimate J5J), and by various rescalings and using 
Littlcwood-Paley theory (as in the work [B], we will discuss this in more detail in 
the sequel) , the estimate (|15|l can be extended to general spherical initial data and 
other L q (L r ) spaces in a straight forward manner. We record this as: 



Theorem 1.3 (Strichartz estimates for spherically symmetric initial data). Let u 
be a spherically symmetric function on M. n+1 such that Du = 0, and set cr^ = n—1, 
then the following estimates hold: 

(16) \\u\\ Lq(Ln < (|| «(0) 11^, + 11^(0)11^0 , 

where ± + 2a < so., and i + a = S _ 7 . 



Remark 1.4. In the language of 0, Theorem 1 1 . 31 savs that the range of indices for 
Strichartz estimates for spherically symmetric initial data are non-sharp (n—1) 
admissible. From the point of view of uniform decay and the general machinery 
of [Hj, this is like saying that these solutions morally decay like t x ~ n . Of course, 
solutions to wave equation in general (if you look everywhere in Minkowski space 
and take the supremum) only decay like t~ %~ . The point here is that if they are 
spherically symmetric, they only do so on a relatively thin set. 
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We would now like to prove a result like Ijl5(l for non-spherical initial data. By 
looking at the Knapp counterexample calculation JJ]J, we see that in order for us to 
avoid a contradiction, we need to replace the term || f[ \\ L 2 on the left hand side of 
that equation by something that is of the order e~ 2 ^~2~T) || f^ || i2 . Now, based 
on our experience with spherical data, it is natural to conjecture that we only need 
to replace the norm || f\ \\ L i with something that incorporates angular regularity. 6 
What we will do then, is test for smoothness of f\ in the angular direction in 
Fourier space. This is done by using the infinitesimal generators of the rotations 
on Euclidean space: 

One sees immediately that: 

(18) ll/illfoi-^ := E H^/ill^ = E W^jffWh ~ WfiWh • 
Interpolating this with the identity, we see that: 

(19) \\tt\w*Li ■= wwnwv ~ e-° , 

where |ri| s = (— A sp /,.)S, and: 

(20) A sph := E n l > 

i<j 

is the Laplacian on the sphere of radius r. 7 Therefore, it is natural to expect that 
if we add s = 2(~ — | — ^) angular derivatives to the left hand side of ©, we may 

get a true estimate in the range 2 ^'_ 2 1 ^ < r. Furthermore, an estimate of this kind 
would be sharp. This is precisely what we will prove in dimensions 4 ^ n, with an 
e loss of angular regularity: 

Theorem 1.5 (Strichartz estimates for angularly regular data). Let 4 ^ n be 

the number of spatial dimensions, and let u be a solution to the homogeneous wave 
equation Uu = 0. Let ctq = n—l be the angular wave admissible Strichartz exponent, 
and let a = be the classical wave admissible Strichartz exponent. Then for every 
< e, there is a C e depending only on e such that the following set of estimates 
hold: 

(21) || u \\ Lq(Lr) < c £ (|| (<ny u (o) + 1| (nyd t u (o) , 

where we have that r + co, s = (l + e )(2=l + |_»=l), I + H = |_ 7; ± + ^;>§ 
, and | + < . All of the implicit constants in the above inequality depend 

6 The regularity in the radial directions does not effect the optimal estimates one gets unless 
you ask for L q (L r ) estimates with q < 2. In order to get these type of estimates, it is necessary to 
incorporate some decay of the initial data, even if it is spherically symmetric. This can be easily 
seen by taking a spherically symmetric wave i\> with the asymptotic , and then forming the wave 

^ = Y*k=o ^(t- k )- Then at time zero one has || ^ (0) || ^2 ~ iV3, but one gets || * \\Li(L r ) ~ N * 
for 2 ^J 2 ^ < r < oo. Therefore it is not possible to take q < 2. 

^Fractional powers of the operator — A Bp h can be defined in the usual way via spectral reso- 
lution. Furthermore, the interpolation identity: (L 2 , — 1 L 2 )t = t L 2 can easily be shown 
using spectral resolution and interpolation of weighted (? sequence spaces. We'll discuss this in 
more detail in just a bit. 



s 



JACOB STERBENZ AND APPENDIX BY IGOR RODNIANSKI 



on n, q, and r. Here: 

\\(nyu\\^ = \\u(o)\\^ + \\\n\ s u(o)\\^ , 

with the analogous norm defined for || (Q) s dtu (0) 

Remark 1.6. A short calculation like the one done above shows that in fact (modulo 
e angular derivatives), all of the estimates 1)21(1 are sharp when tested on Knapp 
counterexamples. Therefore, in this sense, they are all endpoint estimates. It would 
be interesting to try and remove the extra e in these. We will not pursue this issue 
here, although we will do some extra work in the sequel to recover a sharp L 2 
dispersive estimate which could be a start in this direction (see Proposition 15.1(1 . 

Remark 1.7. For the case of n = 2,3 spatial dimensions, we will also prove an 
estimate of the type L 2+ (L°°) and L 2 (L 4+ ) respectively which involves \ an angular 
derivative. However, to obtain the full range of 1(21(1 in the case of n — 3 spatial 
dimensions would require one to prove an L 2 (L°°) Strichartz estimate that involves 
e angular derivatives. While it seems that this type of estimate is out of the reach of 
methods we use here, it should be attainable using the recent method wave packet 
of Wolff ^Hj- In fact, all of the estimates 1(21(1 should be able to be proved directly 
using that method, based on the fact that they correspond to estimates that should 
be true for angularly separated initial data with no extra angular regularity. 

The remainder of this paper is laid out as follows. In the next section, we list 
briefly some of the basic notations we use here. 

In the third section, we list some standard facts about analysis on the sphere 
that will be useful in the sequel, including formulas for Hankel transforms and 
Littlewood-Paley-Stein theory on the sphere. We then use this machinery to re- 
duce the proof of Theorem ll.5l to a suitable set of "endpoint" estimates. 

In the fourth section, we discuss the main tool to be used in this paper: a tp-type 
transform for the Hankel transform. This leads us to consider the localization prop- 
erties in physical space of the bandwidth limited Hankel transform. In particular, 
we provide detailed asymptotics for our Hankel-(/9 transform which will form the 
backbone the Strichartz estimates to be proved here. 

In the fifth, we prove an L 2 dispersive estimate for the wave equation based on 
angular regularity with respect to the momentum operators ((1711 . This can be inter- 
polated with the energy estimate to prove our "endpoints" directly via integration 
in time, avoiding to use of any other analytic machinery such as usual TT* process 
or induction on scales. 

In the sixth, we give a brief description of how our dispersive estimate can be 
modified in a straight forward manner to accommodate multilinear phenomena. 
This shows one strength of the method used here in that one gets the expected 
range of improved multilinear estimates virtually for free out of the machinery de- 
veloped. We will also discuss why these multilinear estimates are sharp, by testing 
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them on a multilinear analog of the Knapp counterexamples introduced above. 

In an appendix to this paper, we provide a proof of the angularly regular end- 
points in the 3 < n regime based on an idea of Igor Rodnianski. The proof there is 
essentially independent of the machinery we develop here (modulo a somewhat sim- 
ilar setup in terms of multilinear estimates) and instead relies on a direct calculation 
involving the energy-momentum tensor for the wave equation. 

2. Basic notation 

We list here some of the basic notation used throughout this paper. For quanti- 
ties A and B, we denote by A < B to mean that A ^ C ■ B for some large constant 
C. The constant C may change from line to line, but will always remain fixed 
for any given instance where this notation appears. We will also use the notation 
A ~ B if there exists a constant C such that ^ • A < B ^ C ■ A. 

For a given function of tow variables, say u(t, x), we denote the mixed Lebesgue 
spaces norms L q {L r ) of u via the formulas: 

\\u\\ q Lq(Lr) := J \\u{t)\\ q Ll dt . 

For a given function of the spatial variable only, we denote its Fourier transform 

as: 

1(0 ■■= J e- 2 ^ x f{x) dx . 

With this normalization of the Fourier variables, the Plancherel theorem becomes 
I / II l 2 = II / IIl 2 j an d one has the Fourier inversion formula: 

f{x) = J e 2 **'* f{£) d£ . 

Using the Fourier transform, we define the homogeneous Sobolev space of order 7 
via the identity || f\\^ := || |£| 7 / ||l 2 - 

For every integer k we define the spatial Littlewood-Paley cutoff operator by the 
formula: 

(22) Kf := Pk f, 

where Pk(0 = Po(2~ h O an d Po 1S a positive smooth bump function, p = 1 on 
the interval (1,2) and zero off the interval (|, 4). With this notation we have the 
following consequence of the Littlewood-Paley theorem for 2 ^ q, r and r < 00: 

(23) IMIL(l-) < E W p ku\\h (Ln • 

k 

For a given function of the spatial variable only, we denote its forward and 
backward wave propagation via the formulas: 

W ± fit,x) = e ±itV ^fix) = J e ±2*m e 2*i*-tf(£) d£ . 
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If u is an arbitrary solution to the homogeneous wave equation, Du = 0, with initial 
data u(0) — f and dtu(0) = g, we may decompose it into forward and backward 
wave propagation in the following way: 



(24) 



u(t, x) 



2i 



J_ (w+g - W-g) + ~ (W+f + W~f) 



Notice that in the case where u = Ui is unit frequency, the functions 

1 1 



hi 



2i 



1 JL 1 f 

7=E 91 + 2 A ' 



2is 



--9i 



have L 2 norm comparable to the H 1 x H 1 1 norm of (/i,<?i). Therefore, in 
the sequel, we will always assume that our unit frequency waves are of the form 



' A hi , and we replace: 
ll«i(0)IU" 



(ll«i(0) + II d tUl (0)||^_0 



3. Some results from analysis on the sphere. 

We list here some basic results from Fourier analysis is spherical coordinates 
which will be used in our proof of Theorem 13.41 We have already introduced the 
two basic differential elements of analysis on the sphere, the rotation vector fields 
{fiij} and the spherical Laplacian A sp h- For every integer $J I, there exists a finite 
dimensional set of functions 3^ on the sphere C M" with satisfy the equation: 

-A sph Y l = l{n + l-2)Y l , Y l e 3>j . 

These sets of functions exhaust the set of eigenfunctions of —A sp h and can in fact 
be identified with the homogeneous polynomials P l on R™ of degree I which satisfy: 

(d?,+^-±d r + ±A sph )P l = 0. 
Setting u> n -i :— |§ n_1 |, and introducing the natural inner product on §™ -2 : 

(F,G) := u-i t f F(uj)H(uj) dtu , 
we have the following basic properties of the Y l : 



Lemma 3.1 (Basic properties of spherical harmonics (see e.g. |10p. 

(1) The dimension of the space y t is \y l \ = \{n + 21 - 2)( n j^ 3 ). 

(2) The spaces yi are mutually orthogonal. That is (Y , Y k ) = for I ^ k. 

(3) Let be any orthonormal basis of y%, then one has the following 
identity for all w G S"" 1 : E» l^/(^)| 2 = 

(4) For each Y l ey h we have the identity: 
Ei<i I! % Yl lli^s-i) =Kn + l- 2)|| Y l \\l 



^From now on, we fix an orthonormal basis {Yj} for each 3^ . For a given function 
F £ L 2 (§" _1 ), we may expand it in the L 2 sense along this basis as follows: 

(25) F = J2^ Y i- 

l,i 



ANGULAR REGULARITY AND STRICHARTZ ESTIMATES 



11 



Using 1)2 5[l. we can define the action of |f2| s on this F as follows: 

(26) p\ s F =^[Kn + /-2)]i^Y/. 

Then using item (2) and (4) of Lemma (|3.1[) . and the fact that 6 3^; (ClijY 1 

is a homogeneous harmonic polynomial of degree I on R n ), we see that we have the 
identity: 

E ll^lli=»(s»-i) = E Z(n + Z-2)|ci| 2 - IIIOI^II^-x) . 

i<j l,i 

Using this, we see that there is, for every test function / on R™ an equivalence of 
norms: 

(27) ll/llini-^ = IIN/IU, , 

where |f2|~ 1 L 2 is the norm from line H18|) . We also make the definition: 

11/ 11^ : = unyfwh , 

:= II /111* + 111^17 Hi* • 

Now, using the fact that Hq is of the form L 2 (l 2 s ) : we have the following standard 
interpolation result (see p^): 



Proposition 3.2 (Interpolation of the angular Sobolev spaces H^). For any set 

or real numbers Si and S2, we have the following interpolation spaces for < t < 1: 

(28) ( #n , ) t = H s n , 
where s = (1 — t)s% + tS2- 

We will also use here some Littlewood-Paley theory in the angular variable 8 . 
We proceed in analogy with (|22|l and let 6 be a smooth bump function such that 
do = 1 on the interval (1,2) and vanishing off the interval (s,4). For each j 6 Z 
we denote its dyadic rescaling as #y(Z) := 9q{2~H). Defining iV = 2 J and using the 
decomposition formula l|25|l . we define the angular frequency dyadic projections of 
a function F on the sphere as: 

(29) F N := E 40 3 (l)Y! . 

We define i<b to be the constant c°, which is the average of F over the sphere. A 
formula similar to l|29() can be used to define /at for a function on the whole of R". 

For a given function F on the sphere, we use line (1) from 13. ll as well as the fact 
that \ (n + 21 - 2)( n +_ i 7 3 ) - when I ~ N to prove Bernstein's inequality for 



"Notice that strictly speaking it will not be necessary for us to use the Littlewood-Paley 
theorem in the angular variable due to our allowed loss of angular regularity. However, the use of 
Littlcwood— Paley cutoffs in the angular variable will be essential for us. 
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the sphere: 

\f n \ < £ \4\-\y!\ , 

l,i : l~N 

< EtE^i 2 ) 1 • N ^ > 

< ( E I4i¥-^ , 

l,i : l~N 

(30) = HFlU^gn-i) . 

By rescaling (|30[) to spheres of various radii, we have the following result on all of 
R n (for N = 0, we of course replace the on the right hand side by 1): 

(31) II IIl'cl-cs?- 1 )) ~ ^^Ik'^/IUw • 

Here denotes the sphere of radius r centered at the origin. 

Next, we record the basic result which allows us to generate certain square func- 
tion expressions in terms of the /n for L r spaces when 2 sC 7': 

Proposition 3.3 (Littlewood-Paley-Stein theorem for the sphere (see ^III |12| . 
and |14|11. Let {0j}°^Q be any set of smooth functions such that there exists a 
< 6 and 6 < 9j(l) < ± for I e (2>',2J +1 ), and 9 3 {l) = for I f {2^\2^ +1 ). 
Furthermore, let 9j satisfy the bounds: 

d M 

|^| < Cm- I- . 

Then one has that for any test function F on S™ -1 , the following bound on the ratio 
of norms holds for 1 < p < oo : 

(32) ^- < || f ii^s-i ) / || (£ |£ ^•(04^r) i IUns- 1 ) < c p , e . 

3=0 l,i 

Here the constant C p ^ depends only on p and the {9j}. 

Using Proposition 13.31 along with the decomposition (|23|l , we have the following 
estimates for functions on space-time, for 2 q, r and r < oo: 

(33) IMIL(l-) < E H^^IIW) ■ 

fcez 

ATe2 N U{0} 



3.1. Reduction of Theorem l3.4l to a frequency localized endpoint. We now 

use the setup we have introduced above to reduce the proof of Theorem 13.41 to the 
following frequency localized "endpoint" estimate: 

Proposition 3.4 (Endpoint unit frequency Strichartz estimate for angularly reg- 
ular data). Let 3 ^ n be the number of spatial dimensions, and let Ui,n be a 
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unit frequency, angular frequency localized solution to the homogeneous wave equa- 
tion Oui^n = 0. Then for every < rj, there exists an < r ri> such that 
r n — > 2 ^Zp as V ~ * and such that the following estimate holds: 

(34) \\ui,n\\l^v) < C v N? +r >\\ui, N (0)\\ L * , 

where the implicit constants in the above inequality depend only on n and r. 

In the case of n — 2 spatial dimensions, we will also prove the following: 

Proposition 3.5 (Endpoint (2 + 1) dimensional unit frequency Strichartz estimate 
for angularly regular data). Let n = 2 be the number of spatial dimensions, and let 
ui.jv be a unit frequency, angular frequency localized solution to the homogeneous 
wave equation □iti.jv = 0. Then for every < rj there is a 2 < q^ such that q v — > 2 
as rj — > and such that the following estimate holds: 

(35) \\u 1>N \\ L «n { L<») < C v Ni + «\\u hN (0)\\ L * . 

Assuming now the validity of Proposition 13.41 we prove Theorem 11.51 as follows: 
Given exponents (q,r) such that ^ + 7 ^ f an d ^ + < we first reduce 
things to the case where q = 2. We define t and r\ via the formulas: 

1 _ t 

q ~ 2 ' 

1 t t-1 



r n I 

Notice that we can find such a ^ i ^ 1 and 2 ^"J 2 1 - ) < 7*1 due to the range of 
(<7, r). Therefore, interpolating with the energy estimate (using Proposition ^. 2|) . it 
suffices to prove <|21ll for indices (2, r\) with 2 ^Zp < r i- 

Next, using the decomposition (|33J) an rescaling the spatial frequency of each 
term in the resulting sum, it suffices to show that: 

(36) || «x,jv Wvpn) < CeN^^"^ \\u hN (0)\\ L 2 . 
To do this, we choose: 

(37) ~ l > - (n - 3) < t e < 1 , 

ri 

such that there exists an < rj with the property that: 

(38) h+V% < (1 + e) 5- • 

2 ri 2 

That such choices are possible follows from our assumptions on the range of r± and 
the identity: 

.1 e. , . ,n — 1 n — 3. 

(5 + 5 *0 = 1+e — , 

2 2 ri 2 

where we have set to — 2 ^"~^ — . Because of the range (|3T)l , we see that is is 

also possible to choose an 2 ^Z^ < r v. ^ 2 ^n-P with the property that: 

t f l-t f 1 



r„. 2(n-l)/(n-3) n 
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Furthermore, using Proposition ^. 41 and possibly a Sobolev embedding, we see that 
we have the estimate: 

2 ( n — 1 ) 

Interpolating this last line with the L 2 (L ™- 3 ) endpoint of we have achieved 
(|36fl . Therefore, in the sequel, we will concentrate on the proof of J33J- 



3.2. The Hankel transform. Finally, to wrap things up for this section, we record 
here the following formula for the action of the inverse Fourier transform on the 
decomposition: 

i.i 

As is well known, this is given by a scries of Hankel transforms. The formula is (see 

my- 

poo ^ 

(39) f{ru) = Y, 2tt(V^I)' / J^l +l {2vrp) c\(p) p^ dp ■ Y*(cj) . 

l,t Jo 

Here J s (y) is the Bessel function of order s. For — j < s, this is given by the integral 
formula: 

< 40 » J -M = m^mmv L e ''" {1 -'^ dt - 



4. The Hankel-i/j transform 

As we see from the formula (|39|l of the last subsection, it is possible to expand 
the expression 9 e -It ^~ A / in terms of spherical harmonics as: 

(41) e-^^/M = Y 4&r)Y!(u) , 

l,i 

where the coefficients c\ are given by the Hankel transform formula: 

POO ^ 

(42) c[(t,r) = 2?r(V^l) 1 / J^ +l (2^rp) e^ 2 ™'" c\{p) p^ dp . 

Jo 2 

Here, as in the previous subsection, the c\{p) are the coefficients in the spherical 
harmonic expansion of f (p). Also, the coefficients c\(t,r) should not be confused 
with the inverse Fourier transform of the c\ (p) . 

We would now like to be able to localize the expressions (|42f> in a manner anal- 
ogous to the localization of the integral (JSJ. This would be a relatively simple 
matter, if we could show that the asymptotic @ held uniformly in n. That is, 
if there was an asymptotic of the form @ for J n-2 . ; (y) which held uniformly as 
( -» co. Unfortunately, it is well known that this is only the case for the region 

^Throughout this subsection, we will work with the operator e _ltv/_A instead of e lt ^~ A . Of 
course this is just a matter of notational convenience, as can be seen for instance by time reversal. 
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I <c y/y (see In the transition regions, that is when ^Jy < I < y, the as- 

ymptotic for J n-2 i ; (y) becomes quite complicated. Roughly speaking, it begins 
to loose oscillations in y while it gains decay in the parameter I. Because of this, 
it does not seem feasible to try and compute an approximate formula for J --2 , ; 
and then substitute it into the integrals l|42[l . Instead we will use a more straight 
forward approach, by first localizing the c\ in frequency as a Fourier series, just as 
we had done for f\ in the integral 10 , and then computing the integral (|42|l directly 
by using appropriate integral representations for the J n -2 | ; . 

Since we are assuming that the initial data in Theorem 13.41 is unit frequency, 
we will assume that all of the coefficient functions c\(p) in the integrals (|42[1 are 
supported on the interval (^, 2). We may take their Fourier series developments on 
the interval (0,4), and we record these as: 

(43) 4(p) = ]T c[ k e^ . 

k 

Expanding the integral l|42|l in terms of the above formula, we see that: 

(44) c\{t,r) = 2n(V=l) l J2 ^ 4,k C|W > 

k 

where: 

/>oo 

(45) = / J^+i V*rp) e~^~> y ( i. 4) (p) dp . 
4 Jo 

In the above formula X/l 4) is a smooth bump function on the interval (^,4). Notice 
that this is not necessarily equal to 1 on any interval because we have absorbed the 
volume element into our definition of 4)- We call the right hand side of H44|) 

the Hankel-ip transform of the function c\ (t, r) . We would now like to be able to 
give a precise bound on how well localized the functions "ipln are m physical space 
for the various values of the half-integer parameter / and the real variable m. This 
brings us to the main result of this subsection: 



Proposition 4.1 (Asymptotics of the functions ip m )- Let ip l m be the function given 
by the formula l|45ll for m = t — j. Then for every set integers ^ N\, N2, there 
exists a constant Cn 1 ,n 2 depending only on the Ni (and the dimension n) such that 
the following asymptotics hold uniform in the parameters I, m, and r: 

(46) 



(47) 

\TpL(r)\ < ? Nl ' N2 , ; ( 7—7 I , 1 < r < |m|+l 

1 mWI (l + r+ |m|)3 • (1 + | r - |m||)^ \ 1 + V 

(48) 

IViWI < — 9 ^r-R{l,m,r) , \m\ + l < r 

(r z — m z ) 2 
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The extra term R(m, r) in line (I48|) above is a positive function with the bound: 

(49) R(l,m,r) < ^ ( (1 + ||m | _ r|)jVl + ^{(7^)^}) ' 

Remark 4.2. The downside of the above asymptotic is of course the region governed 
by Q48JI. When I 2 < r, one can see that the extra factor l|49|) will allow this 
asymptotic to look like (|47|) . Notice that this is consistent with the fact that one 
has the asymptotic @ for J s {y) in this region, and is what is responsible for the 
good localization l|14|l for spherically symmetric waves. Unfortunately, it does not 
seem like one can do much to improve (|48l) in the region where r <C I 2 (except for 
the extra factor of R(l, m, r)). In fact, if one assumes there is an asymptotic for 
V4i( r ) in this region which is of the form (|47J1 . and one sets the c[ k = 1 in the sum 
(I44|l for a fixed l,i, by putting absolute values around the sum l|44(l one would get 
an asymptotic that looks like: |q(0,r)| < -4-. But in this case, c\(t,r) corresponds 
to a delta function along the radial variable in Fourier space, say supported at the 
point Therefore we would have shown a bound like | J s (r)\ < -t- uniform in s. 

This violates the well known asymptotic for Bessel functions: |J s (s)| ~ (see 
|17j). It would be interesting to know if there is a more coherent decomposition of 
the Hankel transform that could eliminate this problem. 

proof of Proposition ^. l\ The asymptotics (|46ll - <|47|) follow more or less directly 
from appropriate integral formulas for the J s {y)- We will need to split the proof 
into the two cases: {r ^ 1 or 1 < r ^ |m| + 1} and {|m| + 1 < r}. 

Case 1: r ^ 1 or 1 < r ^ |m| + 1. 

Here we use a standard integral representation for Bessel functions which differs 
from (|40|l . For s 6 N, one has the following formula: 

(50) J s {y) = I e tycos9 e~ ls9 ale . 

27r Jo 

This can be proved by a simple recursive argument (see |l(Jj . Chapter 4, Lemma 
3.1). All one has to is to show that for both the integral formulas (|50ll and (14011 . 
the function J s {y) satisfies the recursive relation: 

(51) j t [t- s J s (y)] - r s J s+1 ( y ) , o<t, 

for s e N. In light of l|51|l , the equality of (|50() and (|40|l is reduced to showing that 
it is true when s = 0. This can be achieved directly through a change of variables. 

Now using periodicity, integrating over an adjacent interval of length 27r, and 
averaging, we see that for s £ N the following integral representation also holds: 

(52) JJy) = / e tycose e~ tse d0. 

Moreover, by a direct calculation, its not hard to see that the recursive relation 
(|5*T|l is satisfied by both the integrals and (|4*nf whenever s £ \ ■ N. Therefore, 
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throughout the sequel, we may assume that Q52|l is our definition of the Bessel 
function that appears in the integral formula for the function ip l m in dimension 
n. Making this substitution yields: 

/ -\s poo p2tz 



OO J —27T 

2jt 



(53) = 1" x ( i A) (r cos 9- m) e-^+W d9 

4?r J_ 27r u ' ; 



We begin by proving the asymptotic In fact, we will prove a bit more. We 
will show that the asymptotic l|4t)l) holds for r ^ 30. Our first step will be to pick 
up the decay in the (1 + I) parameter by integrating by parts the expression l|53|) 
N2 times. The resulting expression looks like: 

/ • \N 2 N 2 „2tt 

(54) = E / r k Pk {6)x^ A) {rcos9- m )e-^^ e de, 

where the p k (9) in the above formula denote some specific trigonometric polynomi- 
als of degree k who's exact form is not important for our analysis. The next step is 
to gain the decay in r in conjunction with the damping in terms of inverse powers 
of (1 + \m\). To get this, we Taylor expand each jffl A \(h{6) — m), where 

(55) h{9) := rcos9 , 

around the point h = 0. Notice that this is consistent with the fact that we 
are investigating the region where r is bounded. We now define the dimensional 
constant M — p 1 ^], we record this Taylor expansion as: 



2 

M-l 



(56) *S? 4) (ft-m) = E I*£j>(- m ).^ + ± x ^\ u (h) - m) ■ h M 



where u{h) is some smooth function such that |tt(/i)| ^ 30. Substituting the Taylor 
expansion l|56|l into the integral (|54|l . we see that we may write: 

ip l m {r) = A + B , 



where: 



/ ■ \ N 2 , nl-K 

A =(l^- 1 ) E ^X^iH-m) r^ Pk (9)co S i(9)e-^+^d9. 

V — + «/ k+j < M J - 4 ' J - 2 - 



and: 



t , \N 2 N 2 +M 

\B\ < C N2 -r M - hl3T— • E SU P I^^WI • 

We can further estimate the term B above by using the fact that x an d all of its 
derivatives have rapid decay away from the origin. This follows immediately from 
the fact that \ is the inverse Fourier transform of a smooth 0(1) bump function. 
We record this observation as: 



(57) \x)VAy)\ 



( i,4) w! ~ (1 + lyl)^ 1 
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By adding things up and introducing a large enough constant, this allows us to 
write: 

\B\ <£ C Nl . N2 ■ r A ' ■ (^ T] ) • (1 + f M)Nl ■ 

Recalling now that we have set M — f 21 ^] and that we also have 3 ^ n, we see that 
in order to achieve the bound i|46|) . all we need to do is to control the expression 
for A. This is easy to do because a moments inspection shows that in fact one 
has A = 0. This can be readily seen for even dimensions, that is when n is even, 
because in this case the trigonometric polynomials under the integral sign in the 
expression for A are of degree strictly less than ^j^. By orthogonality, the whole 
expression then integrates to zero. In the case of odd dimension, each integral is 
still zero thanks to the fact that is a half integer expression, where as the term: 

p k (9) cos 1 '(B) e~ a , 

is a trigonometric polynomial of integer degree. Since we are integrating over the 
double torus, [— 27r, 2tt], a rescaling turns the expression under the integral sign in 
A into a product of even degree trigonometric polynomials and odd degree trigono- 
metric polynomials. Therefore one has the needed orthogonality. This completes 
the proof of igBJ. 



We now turn our attention to proving the asymptotic 1471) for the regime where 
1 < r ^ \m \ + 1. For the remainder of this section we will assume that to is positive, 
as the other case can be dealt with by a similar argument. Furthermore, using the 
bound that we proved in the previous discussion, we can without loss of generality 
assume that 20 ^ to. 

Our first step is to split the integral on the right hand side of l|53|) smoothly into 
the regions where (1 — cos9) <C 1 and otherwise. To realize this split, we restrict 
the integral (|53|> to the regions: 

Ri = {e\\e\ < i} , 
r 2 = {e\\e±2n\ < i} , 

R 3 = [-2^,27r]\({0||0|<i}U{0||0±2^|<i}) . 

Notice that by symmetry, we only need to consider the regions R\ and R3. On 
these regions, a bit of explicit computation using Taylor expansions shows that: 

(58) h{6) -to = -^r[u(6)} 2 - (m - r) , 6 € R l , 

(59) \h{0) - m\ > y^to , 9 e R 2 and , 
where the function u(9) in the first line above satisfies the bound: 

(60) - < u(9)/9 2 < 2 . 

Notice that to get (|5^|l we have used the condition 20 ^ to. Next, using another 
simple calculation involving the Taylor series of trigonometric functions, as well as 
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the estimates (IKY! . (|5%|) and iffiTijl . we see that one may write for \9\ < 1: 



d9 N ? 



[X{h{6)-m)\ < C N2 sup |x W (^(0)-m)|- E 



^2-^1/11^2-2/0 



(61) 



k=0 



(2 + to - r + r8 2 ) N i+ N i- 



-• E 



A r 2-fe|fl|A r 2-2fc 



1*1 



fc=0 



We are now in a position to bound the integral H53[) on the region i?i . To do this, 
we will employ a smooth cutoff function X(-i,i) such that: 



e (-- -) 



2' 2' 
(-1,1). 



X(-i,i)W =| 
We are trying to prove the estimate 147|) for the expression: 

X(-i,i)(0) ■ Xa A) (h(6) - m) e- l ^+W d6 . 



1 



Integrating by parts -/V 2 times in the above integral and using the bound l|tjl|) (along 
with the fact that we are in dimension 3 ^ n) we compute: 



X(-i,l)(0) ' X iiA) (h(0) - m) e-^+W d9 



1 



l + l 
1 

ITT 



N 2 



N 2 



d N * 

d6 N * 



E 



N 2 



N 2 



X{ -i,i){e)-x aA) {h{6)-m) 
(2 + to - r + r9 2 ) N i+ N i+ 2 



2 1 
f 2 



l + l 



(2 + to - r + 02)iv 1+ Jv 2 +2 

\ A 2 



r N 2 -k^N 2 -2k d0 ^ 

d<9 , 



r 2 (2 + m - r)^ 1 y + 1 



1 + 1*1 



A 2 



(1 + 6> 2 )^ + 

X ^2 



(1 + m + r)2 (2 + m - r) w ± V + 1 / 

Thus, we have proved 147|) for this portion of things. 

It remains to prove the estimate (|47fl for the region R3. We suggestively (c.f. 
(J2U) denote the cutoff here by Xm<\h(e)\ = (1 — X(-i,i))- The calculation is es- 
sentially the same as what was done above, except that here we can afford to be 
more careless about the powers of r which come up through integration by parts. 
Integrating by parts N% times and using the bound l|59|l in conjunction with the 
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estimate (|K7ll . we see that: 

Xm<\ m \(0) ■ X( iA) (h(e) - m) e-^^+O 8 dfl 



< C 



N 2 



l + l 



N 2 



rJV2 ■ , sup , \xfi A) (h-m)\ 



l + l J (i+ TO )*i+ 



_ sup 

A? 2 



iV2 + i 



c 



N 2 



(l + m + r)*(l + m-r) Nl + 
This completes the proof of l|47|l and ends the demonstration of case 1. 



Case 2: \m\ + 1 < r. 

In this section, the main difficulty will be for us to incorporate the remainder term 
(I49[l into our asymptotic. In order to motivate the steps we will take here, we argue 
heuristically as follows: We first define the auxiliary function: 

g(9) := r cos 9 — m . 

It is clear that the main contribution to the integral l|53[) integrating around points 
where the function g{9) defined above vanishes. By symmetry, we need only con- 
sider such points which are positive. We'll call the one closest to zero 9q. Now, 
expanding g(9) around this zero gives: 

(62) h{9) = -<Jr 2 -m 2 (9-9 ) - y (9 - 9 {) f + O((9-9 ) 3 ) , 
where again < 8q = cos -1 ^). Therefore, it is natural to expect that: 

(R.H.S.)H3 « X(i, 4 )(- Vr 2 -m 2 (9-9 )) e~^+ 1 ^ d9 , 



c 



<-2tc 

-i(^+l)6 ,- (n-2)/2 + i „ 

X (kA) {9)e v^CT d6 



(r 2 — m 2 

(r 2 - m 2 )2 V Vr 2 -m 2 J J-oo u ' 1 
where the X/Z^ > ^2 = 0, ±1, ±2, . . ., are the (unique) set of derivatives and deriva- 

(.4 I 4 ) 

fives of X(i 4) which vanish at infinity. These satisfy the bounds similar to (|57|1 . 
even if JVa is positive: 

(64) x^/My) < 



HiA) yy > ^ (l + \y\) 2 ■ 
This, of course, comes from the fact that X(i 4) is the inverse Fourier transform of a 
unit frequency function, and is therefore supported away from the origin in Fourier 
space. Using (|64|) we easily get the following bound for arbitrary non- negative A^: 



|(r.h.s.)E3| < ( Aj 



(r 2 — m 2 ) 9 




±N 2 ' 
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In order to make the previous argument rigorous, we need to justify the two 
approximations used in the lines directly above (|63|l . We will not be able to do 
this completely, which will be responsible for the extra term in the formula 149(1 for 
R(l,m,r). As in the previous two subsections, we argue by isolating the interval 
of integration. By symmetry, and the decay bound 1|57JI. we can without loss of 
generality assume that we are integrating (|53[) over the interval [0, ^p]. We will 
first need to go a little further and chop some more off of the left hand side of this 
interval. What we'll do is take 6\ (we're still keeping 9 — cos~ 1 (^)) to be such 
that: 



(65) 



: sin 







Taylor expanding sin 9 gives us the bound: 

(r — m)i 



(66) 



< 



7" 2 



Furthermore, we also have the bound: 



(67) -(r-m) ^ 



-m < rcosO-m = g(8) , 6 <E [O,0i] 



Now, using i|66|) and (|67|l above, in conjunction with the asymptotic l|57|) . we im- 
mediately see that: 



|X(J,4)G?(0))| d9 



< 



(r 



C 



r i {r-m) N ^ ' 

This is enough to give l|48|l for this portion of things. Therefore, it remains to 
compute the integral: 



h 



X(i,4)(#(0)) 



'-+1)9 



(10 



Keeping the Taylor expansion l|62l) in mind, we now make the following change of 
variable for 7i: 

1 



<p{B) 



VV2 



Using this, we can write: 



(68) 



I'll = 



1 



V7 1 



*(i,4)(¥»(0)) 



; (n-2)/2+t 



(16 



Where R is the interval: 



R 



y/r 2 - m 2 (9 1 - 6 ), \f7 2 



— 2/ 3?r 
4 



%) 



The desired result will now follow by integrating by parts as many times as necessary 
the integral l|68() . However, some care needs to be taken in order to control terms 



involving <p' , which can be as big as 



Also, one needs to know that the 



\J r 1 —m,' 2 

higher derivatives of ip possess some decay in order to deal with terms where the 
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derivatives fall on ^7 instead of the exponential factor. We address these issues 
now. The first observation we make is that we may integrate the bound: 

(69) \ < / 2 sinfl, e [9x,^] , 
using the fact that v?(0) = to get that: 

(70) \e\ < \<p{9)\ . 

over the range of integration. However, this is not enough to address the issue of a 



product of the form (tp') k ■ xfl *y To handle this, notice that one has the bound: 



(71) r sinfl < m-r cos6 + l, € [61, ^] , 
V r 2 — m 2 4 

as can be seen from the fact that this bound is true for 6 = 0i, and that upon 
differentiating both sides of this expression one is reduced to showing (the inequality 
for the derivatives is trivial for S [§, ^f- ]): 

1 7T 

— === < tan0 , e [0i, -] . 

\/r 2 — m 2 2 

This then follows immediately from the increasing nature of tan on this interval 
and from the identity tan(0i) = ^ r ' 2 ~ m2 , Now, combining the bounds (|71() . I|70|) 
and l|57|). we see that for any positive integers k,j, one has the estimates: 

(72) •*<?>(«)) ^ (ifjlrp- 

Finally, we record the fact that the higher derivatives of y satisfy the following 
simple bounds, which can be verified through direct calculation: 

(73) M fe) | < 1 —j- T , M = 2,3,... . 

(r — m) K 



We are now ready to bound l|t)8|) . First, we integrate as many times as necessary, 
letting the derivatives fall on the term ^(i 4)(y(0))- The resulting expression looks 
like: 



< 



Vr 2 



X(i,4)(v(#)) e 

— \ N 2 



; (»-2)/2+i 



(20 



1 + / 



d^ 2 



f/0 



E i^s*(i.4)(*»w)i 



(77? 



Expanding out the derivatives in the first term of the above expression, and using 
the bounds l|71|l and (|70|l , we see that we can bound: 



^X(i,4) (*>(*)) 



< 



Cn 2 



(l + l^l) 2 



Therefore, the integral in the first term above does not cause us any trouble. As 
far as the boundary values are concerned, we can use the bounds l|65fl. IJB7J, (|7T|l . 
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and |K7| to show that: 

E 1^(4^(^)1 



fc=0 



OR 



<C N2 sup (b (fel) W| fe2 -|xfi)(^)) 



v / r 2_ m 2 ( 3 f _ So) 



v / r 2 -m 2 (6»i-e ) 



< c Nll 



N 2 



1 



1 



(r — m)-^ 1 (r + to) 



< 



Cai, 
(r - to)^ 1 



N 2 



It remains to bound l|t)8l) where we integrate by parts and let the derivatives fall 
on the exponential factor. Doing this N2 times yields: 



X(I,4)(<^)) 



.- (r.-2)/2 + i 

e 



d0 



iV 2 



< 



(74) 



ip' (9) d9 



(42 ) (^)))- e " i ^ 



(n-2)/2 + i 



d0 



1 



N 2 



E i*lw*>) 



</(6>) d0 
, 1 



JV 2 



■- ("-21 



99' (9) \<p'(0) d9 



fc-i 



d0 



,(n-2)/2+i 



e 



To control the first term above, we bound: 

A 2 



(75) 



f 



< 



<p'(9) d9 



sup 

0<Ck<2N 2 -l 



■ (n-2)/2 + i 

e v^CT 
JV 2 



1 _ 

—\ k -Y 



1+1 



\/r 2 — to 5 



|iV a -i(v (1+,) W)| , 



where Po is a constant, and the other denote a homogeneous expression of 
weight k in the variables (</^ 2 \ . . . , ip^\ ■ . .), where each y>W is given weight i — 1. 
Therefore, using the bound (|69|l (which in particular implies that |1/V'| ^ 2) as 
well as the bounds (|73|l . we see that we can estimate: 



N 2 



CZ3 ^ Cjv 2 E 
3=0 



1 + / 



\Jr 2 — m? 



1 



< CV 2 



1 + / 



iV 2 



(r - m) w 2-3 
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Finally, to deal with the boundary terms on the right hand side of Ij74(l we use the 
bounds (|H7| and and simply estimate: 



N-2 



k=l 



if' (6») W(0) dd 



k-1 



„. (n-2)/2 + i 

e 



<C N , sup |*5i*A(^))| 

0^k^N 2 U ' ' 
CjVi ,N 2 



OR 



< 



OR 



Combining the above estimates together, we have shown that: 



Hi I < 



1 



Cn u n 2 ( 

Vr 2 - m 2 \ (r- m)^ 



mm ■ 



I 



±N 2 ' 



This completes the proof of the asymptotic l|48(l , and the demonstration of Propo- 
sition O □ 



5. Some L 2 dispersive estimates for the wave equation; linear and 

BILINEAR ESTIMATES 

We are now ready to directly proceed with the proof of estimate Ij.'Mfl . As we 
have mentioned previously, we may assume that mi.jv is of the form e -**v /3 ^ /i.jv, 
for some unit frequency and dyadic angular frequency function /i,j\r- For this iti.jv, 
we the formulas from lines (|4 1 11 - l|42(l . H43|l . (|44|l . and (|45|) to expand its wave 
propagation into harmonics and Hankel-t/? transforms as follows: 

(76) e' u ^h, N (r) = £ c{ k ^_,(r) • Y* 

l~N 

In the above expression, we have absorbed the constants 27r(v / — 1)' into each c\ k . 
By orthogonality of everything in sight, we have that: 

(77) ~ E i c Ui 2 - 

i,k 
l~N 

We now state our main result as follows: 

Proposition 5.1 (1? dispersive estimate for the wave equation). For dimensions 
2 n, let e-^v^ j be 

given by the formula (|7t)|) . Then one has the following 
estimate uniform in t and r: 
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Before proceeding with the proof, let us first show briefly how the above estimate 
may be used to show Ij34|l - H35|l . We begin with (jj2J. What we need to do is to 
show that for every < rj, there exists an 2 ^Z^ < r v> sucn that the following 
estimate holds: 



(79) 



-it\/-A 



h,N\\mL^ } < N^-\\f hN \\ 



where the implicit constants depend on both 77 and r. We also need r n to approach 
2 ^n-2 ^ as V ~ * ^' Now, interpolating 10 (f TH|> with the energy estimate (j2J gives the 
following result for 2 ^ r„: 



(80) 



-iW-A 



E 



J |2 



■2V V 



Choosing r, ( according to the formula \ + i] — ( | — i) 2^1 , we see that we have 
the following identity holds: (| — ^-)2(n — 1) = 1 + 277. Therefore, for this choice 
of r v , we may square (|8UII and integrate directly in time to achieve (|79|l . Also note 



that r„ — > 2( -" ^ as r? — > 0. 

'/ n— 2 * 



In the case of estimate (J2HJ, that is n = 2 spatial dimensions, we use Holders 
inequality to see that for every < 77 the following estimate holds: 



< 



|2 



i.k 



;e i c Ui 2 )^ 



Therefore, of we choose = 2(1 + 77) a direct integration and the energy bound 
(|77|l shows that we have: 

\\e- U ^* fl <N \\ L «n(L°<>) < ■ II fi, N \\ L 2 . 

We are now reduced to proving the estimate (|78|) . 



proof of Provosition \5. 11 Writing everything out, and using a Cauchy-Schwartz in 
the sum over i for each fixed / we have that: 



(81) 



h,N Ml 



l~N 
i.k 



< 



ElE r ^ • su p E \ y i 



i l~N 
k 



l~N 



Ik? 



10 This can be achieved by simply interpolating in weighted i 2 and the usual Lebesgue spaces 
using the map {c\ k } — > e _It%/-A fx ff for fixed time (see 0). 
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Now, using the formula in item (3) of Lemma 13.11 and the fact that j (n + 21 
2 )( n tli 3 ) ~ N "~ 2 for I ~ N, we have: 







(L.H.S.)JHU ^ 


EE 




\ k 




( 


< 


EE 




^l~N k 



r 2 



r ~ 2 chipLdr) 



TV- 



TV— 



where the last line above follows from a Cauchy-Schwartz and the bound I 
The proposition will now be shown if we can prove that: 

J 12 

(82) x * — 1 ' 1 1 ' 1 ' x * ' l ' kl 



N. 



E^ ^ (E 7i 



(i + i* -ID"- 



i.k 



We now use the asymptotics of Proposition ^. II and a little computation to split: 
where b(t, t — |^ is supported where 1 < \t — || < (r — 1), and: 



A[r,t- 



(i + i^ii)^ 1 (i+k-|t-|i|) ' 



(i + lt-ll)^ (i + k-|t-fi|) 



T -R(l,t- ~,r) . 

t- 4 hs V 4' / 



Here R,(l,t— \ ,r \ is the function from line l|49l) . Substituting A^r, t — ^ into the 
left hand side of the sum 182|) . and using a Cauchy-Schwartz immediately yields 
the desired result for this half of things. In order to finish up, then, we only need 
substitute to B-(r,t— jj into the right hand side. Doing this and using a Cauchy- 
Schwartz and re-indexing, we see that it is enough to show: 

<w> E d + i-'-im) - *('■*■" si- 

fcGfZ+t V 1 ' l|y 
|fc|<(r-l) 

uniform in r and I. This can be done by comparing things to the appropriate 
integrals and is left to the reader. Notice that the convergence factor R(l,k,r) 
avoids any logarithmic divergences in this sum. This completes the proof of i|78|) 
and therefore the proof of Theorem 1 1.51 □ 



6. Bilinear estimates for angularly regular data 

As is well known, for applications in the lower dimensional setting, linear esti- 
mates of the form (|21[) are not sufficient. What is needed are multilinear versions. 
An extremely versatile method for building these type of estimates is based on the 
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"fine and coarse scale" machine of T. Tao (see 0). The basic idea is to fix a scale, 
say -jj for /j < 1, and then decompose the domain of spatial variable into cubes 
with side lengths of this scale. Then, one replaces the usual L r norm in the spatial 
variable with £ r (L 2 ), where the L 2 norm is taken on the "fine" scale of each indi- 
vidual cube, while the £ r norm represents the "coarse" scale which is summation 
over all cubes. One reason this method is so important, is that it allows one to use 
the bilinear construction process directly in an iteration procedure where resorting 
to the canned estimates that this method ultimately provides may be unduly bur- 
densome. Also, the way these estimates are constructed will allow us to generate 
bilinear estimates where only one term in the product contains extra angular regu- 
larity. These type of estimates are extremely important for applications when one 
proves inductive estimates via an iteration procedure where it is necessary to con- 
sider estimates for products with an angular derivative falling on one or the other 
term in the product. Because of these considerations, we will content ourselves 
here with the dual scale estimates themselves, and not bother with listing out the 
various multilinear estimates which follow from them. In the (4 + 1) and higher 
dimensions, the usual dual scale estimates read: 

Theorem 6.1 ("Improved" Strichartz estimates). Let 4 ^ n be the number of 
spatial dimensions. Let < fi < 1 be given, and let {Q a } be a partition o/K™ into 
cubes of side length ~ — . Then if u% is a unit frequency solution to the equation 
□ill = 0, the following estimates hold: 

(84) IKEII^WII^wJ'N ~ A* -1 II «x(0) ]U- , 

a 

where ^ < r. 
n— o 

We show here that both the range of (|84(l and the power of \i that appearers there 
can be significantly improved. In what follows, we will only bother proving an es- 
timate which is the analog of (j3.4|l . More general estimates which involve various 
amounts of angular regularity can then be gained by interpolating this estimate 
with (|84|l . Also, it is not so easy to work out the Littlewood-Paley theory in the 
angular variable for localized norms like those that appear on the left hand side of 
(I84|l . However, since we are already loosing a small amount of angular regularity, 
one may simply replace the Littlewood-Paley sum lj33J) by an i 1 sum. As we have 
noted before, this has no bearing on applications. 

Theorem 6.2 ("Improved" frequency localized Strichartz estimates for angularly 
regular data; endpoint case). Let 3 ^ n be the number of spatial dimensions, and let 
Ui,jv be a unit frequency and angular frequency localized solution to the homogeneous 
wave equation □wi 1 at = 0. Let < /i < 1 be given, and let {Q a } be a partition o/M" 
into cubes of side length ~ — . Then for every < 7], there is a and 2 ^_P < r r) 

depending on r\, such that r v — > ^iB— U as n — > such that the following estimate 
holds: 

(85) \\ C£\\ ui , N (t) \\%( Q J^ hi < C^-i-^N^WuxAVh* ■ 
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Remark 6.3. Up to the small loss in ^ and N, the estimate (|85|) is sharp when 
tested against the bilinear analog of the Knapp counterexamples (J5J. We construct 
these briefly as follows. We consider (frequency) initial data sets \B' which along 
with being highly localized in the angular variable, are also well localized in the 
radial variable. That is, we are now assuming \B' is supported on a small square 
of dimensions ~ e x e x . . . x e, lying along the £1 axis between 1/2 < £i < 2. A 
quick calculation then shows that the integral: 

behaves like \e lt ^~^f{ (x)\ ~ e n on the space-time region Sl x : 

t = 0(e- 2 ) , 
t + xx = Oie- 1 ) , 
x' = O(e^) . 

Choosing our cubes Q a with side lengths — ~ i, we see that for any 2 ^Zp < r ( m 
fact for any 1 ^ r ^ oo) the following is true: 

(86) f,-ie-i\\fl\\ L , ~ e"-i-i < || || e «^/J 11^))" U 

Q 

Therefore, using l|19|) . we see that (|85|l is indeed sharp for this sequence of initial 
data. Of course the condition ^BrlJ < r cannot be improved, even for spherically 
symmetric initial data. In fact, one can also see for these type of waves (say with 
the asymptotic (7J), l|85|) is again sharp. We leave this simple calculation to the 
interested reader. 



proof of estimate l|85|l . The proof will be similar in spirit to that of 1)34(1 ■ However, 
we will need to use orthogonality in a more fundamental way here. This is because 
we will not be able to rely solely on L°° as we did in the proof of (|78|l . Ultimately, 
this has to do with the fact that the wave packets ip l m (r) are not well localized in 
physical space when \m\ < r and r<I 2 (see the remark after Proposition 14. lj l. 

In order to prove estimate Q85[). we begin by writing the basic energy estimate 
(J21 in the following way (we are still using the notation from line (JZBJ): 



(87) 



-iW—A f 1 1 2 

h,N\\ 



< 



E 

, i,k 



U 2 



Our next step is to prove the following spatially localized fixed time estimate: 

/ 



i-l n—l 



sup \\e-^-*f hN \\ L 2 {Qa] < \\n^- — N— 



•sr^ \ c i,k\ 2 

h (i+i*- IP"- 1 
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Assuming for the moment the validity of Ij88(l . we may interpolate between H87JI 
and l|88|) to achieve the following estimate for 2 ^ r v : 



(89) £ \\e-^A, N \\Z iQo 



< 



E l H I Tn D I -Un^-^iV 



-l 1 1 l 



.k 



As in the previous section, choosing r v by the identity, ^+f] = (^ — ^r) ^^"j we have 
that (i — — )2(n — 1) = 1 + 2r]. Therefore, squaring (|59"|) and integrating directly 
in time, we will have achieved (|85() . Therefore, we now concentrate on proving (|88|l . 

It suffices to show (|88f) for a fixed a. Therefore, we will now assume that we are 
on a fixed cube Q a . Our first step is foliate Q a with the hypersurfaces D Q a , 
where §™ _1 is the sphere of radius r centered at the origin. Notice that for each 
r, one has the area estimate | D Q a \ < Therefore, using a Cauchy- 

Schwartz, it suffices to prove the estimate: 
(90) 



^ UV ^h,N\\ Ll{L ^-. )){Qa) < |ln M |iV^ \J2 7j 



Id J 2 



; .V (1 + 1*- ID"" 1 
i,k 



Next, we chop Q a into at most | ln/i| dyadic pieces, which are of the form Rj D Q a , 
where Rj = {x | 2 J < |a;| < 2 J+1 } is the radial dyadic region of size r ~ 2 3 . We 
only need to do this for < j, i.e we keep the ball of bounded radius (say radius 
= 2) as a single region R . Therefore, to show (|90|) . it suffices to estimate: 



( 91 ) »P lle-^/i^ll^L-W-ww S ^ E (i + it-lijn-i 

Now fix Rj. To prove estimate (I91f) here, we run the decomposition l|7()l) and 
decompose the resulting sum into the sum of pieces: 



-itV-A 



/i,iv — S + Si 



where: 



i,k 



si = E r ^ <*v»U( r )- y < 



|t-*|^4-2 J 
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To estimate on the sum Ei, we simply keep it as a whole object and use the 
Bernstein inequality to estimate: 

Si 



sup 

i 



< N— 2" 



Il2(loo(§: 

^IISilU, 



< N— 2 — ~ 



E 



J |2 



< N~ 



E 

l~N 



(i + l* 



II)"- 



/ 



As was to be shown. Therefore, we are reduced to bounding Eo- Notice that for 
each fixed r £ Rj, and for each term in this sum, we have \t — j\ > r. Therefore, 
we can use the well localized asymptotic (|47f) and a Cauchy-Schwartz in the k 
summation (in a way similar to the computation started on line Q81|l) to bound this 
term in L°°(8™ -1 )) by: 



Mr) 



< 



( 



V 1 



E 

i,k 
-|I>4-2 J 



(1 + 1* -II) 7 



(1+ r-|t-|| ¥ 



N- 



J 



Squaring this last expression, and integrating each term in the sum with respect to 
dr (notice there is no extra volume element because we took the sup on we 
see that we have proved (|91|l for this portion of things. This then completes the 
proof of estimate (|88p. and therefore the proof of l|85|) . □ 



7. Appendix 

We present here a simplified proof of Theorem 13.41 and Theorem 16.21 First, we 
note that in order to prove both (|31Jl 

(92) \\UxM\LlLl^ C V Nh+V \\ U hN{G)\W 

and JSSJ 

(93) \\uiM\Ll^Li<Cr l ^- r, N 1 - +1 'i\\u 1 , N {0)\\ L 2 

with p v — > 2(n — 1)/(n — 2 — rf) and arbitrary small 77 > 0, it suffices to prove the 
estimate (I85f) . Here ui.at is a unit frequency solution of the wave equation Du = 
of angular frequency N and the norm 
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where {Qa} is a partition of K n into cubes Q a of side length of /it -1 . 

To see this, we choose a partition {Qa} of size fj, = 1 and compute, using the 
Sobolev embedding on Q a , that for any 2 < p: 



II Ui, N \\ L?L p = || Yjui.jv 11 '' 



LP(Q a ) 



< ii EiKwiis 



L 2 (QJ I Nit ■ 

Therefore, we see that it suffices to deal with the estimate IjHHJI ■ First, using the 
Sobolev inequality C (fl}^ n ^~ if, on the unit sphere S ra_1 for angular frequency 
localized functions, we have the following estimate for any tiling of K™ by cubes 
{Qa} of side length ~ 

IK,iv(t) ll£ 2 (<2«) ~ / Wi^it^dx^fi-^ drr 11 - 1 sup K,Ar(rw, t)| 2 

Jruj£Q a JO uES"- 1 

< || (1 + r)-^ (Q)^ Ul , N (t) \\ Ll . 

Thus, 

£ fi-^N 2 ^ 1 || (1 + r)" 2 ^ ui jAr (t) |U|. 
Interpolating this with the trivial estimate 

||ui,w(i)||£ 2 ij < || Ui tN (t) \\ L 2, 
we arrive at the following estimate which will be our point of departure: 

(94) < C„ ^-*- 7J || (l + r)-i( 1 +"><0i( 1+ ^)w l! ^|U ?jLS , 

where p v = ^-2^2 ■ Therefore, using and that (f2) 3( 1+? ?) Ul!Ar is another unit 
frequency solution of the wave equation localized at the angular frequency N, we see 
that in order to prove (|93[) . it suffices to prove the following space-time Morawetz 
type estimate for unit frequency solutions of the homogeneous wave equation when 
3 ^ n: 

(95) ^ [ 1 \ Ul (t,x)\ 2 dxdt < ||m(0)||| 2 . 

In the above estimate, the implicit constant depends on rj and the dimension, and 
from now on we will keep this dependence implicit. Estimate l|95l) . as well as as 
its version with r\ = and a logarithmic loss in time, was proved in the work 
of Keel-Smith-Sogge [5] in dimension three with the help of the sharp Huygen's 
principle. We will prove (|95|l directly using essentially nothing but an integration 
by part argument. Before we continue, let us make one more reduction. It turns 
out that l|95|l is more naturally proved for the spatial gradient V^ui. Since u\ is 
unit frequency, this reduction does not effect the validity of l|95|l . Indeed, suppose 
we have the following estimate: 



(96) r f 1 \V xUl (t,x)\ 2 dxdt < ||V MUl (0)|| 2 L2 
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for any unit frequency solution u\. Let u\ be another unit frequency solution with 
the property that Aui = — u\. Then for any e > 0, 

/ / t- rrr - \ u i(t, x)\ 2 dxdt = I / — -r— A{*iui dxdt 

f°° f 1 f°° f rV r 

= J ^ ^ ( TT7 ^V^ 1 V^ 1 dxdt + (! + ,) J ^ -—^V,** dxdt 

1 f°° f 1 

'V x ui(t,x)\ 2 dxdt + / — — ■ — -7— |V x wi(t) x)\ 2 dxdt 



f 



and l|95l) follows from (|96[) by choosing a sufficiently small e. Therefore, we may 
now assume that we are trying to achieve estimate l|96[) . The fact that H96|) con- 
tains the gradient V ' x will allow us to prove it for an arbitrary solution to the 
homogeneous wave equation. This will be done using some more or less standard 
energy-momentum tensor techniques. Let 4> be a solution to the homogeneous wave 
equation: 

(97) a<f> = , 
and let Q a p[4>\ be its energy-momentum tensor: 

QapWi = d a cf)dp<f> - -g a /3d y (j)d^(p . 

Here the Greek indices run on the set a, /3,'y = 0, . . . , n, and g a p = diag(— 1, f , . . . , 1) 
is the standard Minkowski metric. The key feature of Q a p[<fi\ is that it is space-time 
divergence free: 

(98) D a Q aP [4>] = , 

where V denotes the Levi-Civita connection of g a p. We now contract Q a p[4>] a 
radial vector-field: 

X = f{r)d r . 

Define the momentum density: 

P a [<f>,X\ = Q a[} [4>]xP . 
A quick calculation shows that the divergence of P a [<j>, X] satisfies the identity: 

(99) D a P a [^X] = l -Q a ^]^ , 
where: 

TT a f3 — DaXp + DpX a , 
is the deformation tensor of X. Introducing the orthonormal frame {dt, d ri ca}, 
where the {e^} form an orthonormal frame on each n— 1 sphere S^ r = {?' = const.}, 
one can easily calculate this quantity to be: 

/0 ... 0\ 

(100) 7r 0/J = f'(r) ... 

\0 ~r~$AB 
where Sab denotes the metric on the spheres St. r - In particular, we see that: 

(101) to = /'(r) + (n-l)^ , 
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Next, using (|99|l and 1(100(1 . a direct computation shows that: 

2D a P a [<j ) ,X] = /'( r )|a r 0| 2 + ^^| 2 -itr7r9^a 



In the above formula, ^(f>\ 2 = 8 AB \eA{4>)\ ' |e s (0) | denotes the angular portion of 
the spatial gradient |V X </>| 2 - Since, d 1 4>d 1 4> = ^l^l 2 ) if we define the modified 
momentum density: 

?*[</>, X] = P a [<f>,X] + jtnr<t>d a cl>-^d a (trir)\<f>\ 2 , 
we end up with the identity: 



(102) D a P a [^X] = - ^f>( r )\drW+ J -j^W<l>r) - gA(tar) \<t>\ 

Integrating ((102(1 over a time slab, we arrive at the following a-priori estimate for 



(103) / P [(/>,X](Q)dx = / P [ct>,X]{T) dx 



\ (f'(m4\ 2 + ^\f<t>?) - (tnr) |0| 2 dxdt 



dx 



where: 

(104) / P [<P,X}(0) dx = f (x a a t 0(O)a a 0(O) + itr7r0(O)5 t 0(O) 

with an identical expression for the time = T boundary piece on the right hand 
side of 1103|) above. In particular, using (|101fl . one has that if \f(r)\ < 1 then: 



P [(j),X}(0) dx 



(105) < ||V t ,^(0)||i 3 , 

where in the last line we used the Hardy inequality. Similar estimate holds for 
the other boundary term P [</>, X] (T) . Therefore, using the usual conservation of 
energy, we have the following a-priori estimate for solutions to the homogeneous 
wave equation: 



(106) | J J - ^f'(r)\d r <P\ 2 + ^W<P\ 2 ) - gA(trTr) |0| 2 dxdt | 

< ||v t ,^(o)||| s . 

We now use (jlQ6fl to derive (|96fl by choosing the weight function f(r). 

r 

f(r) = — - 
e + r 

for some e > 0. A direct calculation shows that, with this choice of /, we have 

. , . / e n — Is 1 f , , . . er 3e 2 (n — 1)\ 

A (tor = A( r — — ) = - ( (n-3 r+3 n-3 — '-) < 

(e + r) z e + r r{e + r)°\ e + r e + r / 

in dimensions n > 3. Therefore this gives the a-priori estimate: 

(107) r [ (-^-^\dA\ 2 + ^-^^\ 2 ) dxdt < ||v M 0(o)|| 

.l—^.lvn \{e + r) e + r / 
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with an implicit constant in < independent of e. In particular, choosing e = 2 k 
with an integer k > 0, we obtain 

/ / IV^I 2 dxdt+ [ [ dxdt < ||V tjX 0(O)||£ a . 

J~ooJ\x\<l J-oo J2 k - 1 <\x\<2 k + 1 r 

Dividing the above inequality by 2~ 1) ' c and summing over k immediately yields the 
desired estimate 
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